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Abstract 

The framework of relativistic self-consistent mean-field models is extended to include correla- 
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. tions related to the restoration of broken symmetries and to fluctuations of collective variables 

The generator coordinate method is used to perform configuration mixing of angular-momentum 
' and particle-number projected relativistic wave functions. The geometry is restricted to axially 



symmetric shapes, and the intrinsic wave functions are generated from the solutions of the rel- 



O \ ativistic mean-field + Lipkin-Nogami BCS equations, with a constraint on the mass quadrupole 

^3 



moment. The model employs a relativistic point-coupling (contact) nucleon-nucleon effective inter- 
action in the particle-hole channel, and a density-independent ^-interaction in the pairing channel. 



> 

X. 

, Illustrative calculations are performed for ^^Mg, ^^S and ^^Ar, and compared with results obtained 

employing the model developed in the first part of this work, i.e. without particle-number pro- 
jection, as well as with the corresponding non-relativistic models based on Skyrme and Gogny 
effective interactions. 

PACS numbers: 21.60.Jz, 21.10.Pc, 21.10.Re, 21.30.Fe 
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I. INTRODUCTION 







In the first part of this work [l] we have extended the theoretical framework of relativistic 
self-consistent mean-field models to include correlations related to the restoration of broken 
symmetries and to fluctuations of collective coordinates. In the speciflc model which has 
been developed in the generator coordinate method (GCM) is employed to perform con- 
figuration mixing calculations of angular momentum projected wave functions, calculated in 
a relativistic point-coupling model. The geometry is restricted to axially symmetric shapes, 
and the mass quadrupole moment is used as the generating coordinate. The intrinsic wave 
functions are generated from the solutions of the constrained relativistic mean-field + BCS 
equations in an axially deformed oscillator basis. In order to test the implementation of 
the GCM and angular momentum projection, a number of illustrative calculations were per- 
formed for the nuclei ^^^Hg and ^^Mg, in comparison with results obtained in non-relativistic 
models based on Skyrme and Gogny effective interactions. 

In this work we develop the model further by including the restoration of particle number 
in the wave functions of GCM states, i.e. we restore a symmetry which is broken on the 
mean-field level by the treatment of pairing correlations either in the BCS approximation, or 
in the Hartree-Fock-Bogoliubov (HFB) framework. We perform a GCM configuration mixing 
of angular-momentum and particle-number projected relativistic wave functions. Projection 
on particle number is crucial whenever the number of correlated pairs becomes small and 
the density of levels close to the Fermi energy is low, a situation typical for the description 
of phenomena related to the evolution of shell structure : reduction of spherical shell 
gaps and modifications of magic numbers in nuclei far from stability, occurrence of islands 
of inversion and coexistence of shapes with different deformations, moments of inertia of 
superdeformed bands, etc. We thus plan to build a self-consistent relativistic mean-field 
model in which rotational symmetry and particle-number are restored, and fiuctuations of 
the quadrupole deformation are explicitly taken into account. Such a model can be applied 
in a quantitative description of shell evolution, and particularly in the treatment of shape 
coexistence phenomena in nuclei with soft potential energy surfaces. 

In Sec. [TTl we outline the relativistic point-coupling model and the Lipkin-Nogami ap- 
proximate particle number projection, which are used to generate the intrinsic mean-field 
wave functions with axial symmetry, and introduce the formalism of configuration mixing 
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of angular-momentum and particle-number projected wave functions. In Sec. Illll the model 
is tested in a detailed analysis of the spectra of ^^Mg, ^^S and ^^Ar. In order to illustrate 
the effects of particle-number projection, we compare the results with those obtained em- 
ploying the model developed in the first part of this work l]], in which the intrinsic wave 
functions are generated from the solutions of the constrained relativistic mean-field + BCS 
equations, without particle-number projection. The results are also discussed in comparison 
with the corresponding nonrelativistic GCM models based on Skyrme and Gogny effective 
interaction. A brief summary and an outlook for future studies are included in Sec. [IVl 

II. THEORETICAL FRAMEWORK 

A. Implementation of the Lipkin-Nogami pairing scheme 



In the model that we have developed in Ref. the intrinsic wave functions are generated 
from constrained self-consistent solutions of the relativistic mean-field (RMF) equations for 
the point- coupling (PC) Lagrangian of Ref. ^. Only basic features of the RMF-PC model 
are outlined in ^|, and we refer the reader to and references therein, for a complete 
discussion of the framework of relativistic point-coupling nuclear models. The specific choice 



J] defines the mean-field energy of a nuclear system 



k 

+ 



of the PC Lagrangian [l, 4 

Ermf = y £RMF{r) 

^ dr vl tpkir) {-i-fV + m) ipkir) 

( + yP5 + -jPs + yP5Ap5 + + y OmJ'^) + yJ^Aj'^ 

+ -^jTvijTv)^. + —J^v^Utv),. + -I^PtS + —PTS^PTS + l^PpA 1 , (1) 

where "0 denotes the Dirac spinor field of a nucleon, and the the local isoscalar and isovector 
densities and currents 

Ps{r) = tpk{r)Mf) , (2) 

k 

pTs{r) = ipk{r)T3ip^{r) , (3) 

k 

f[r) = J2^lMr)rMr) , (4) 



JTvir) = Y.'^l Mr)rr,^lj,{r) , (5) 

k 

are calculated in the no-sea approximation: the summation runs over all occupied states 
in the Fermi sea, i.e. only occupied single-nucleon states with positive energy explicitly 
contribute to the nucleon self-energies, vl denotes the occupation factors of single-nucleon 
states. In Eq. ([T]) pp is the proton density, and denotes the Coulomb potential. 

In addition to the self-consistent mean-field potential, for open-shell nuclei pairing cor- 
relations have to be included in the energy functional. In this work we do not consider 
nuclear systems very far from the valley of /5-stability, and therefore a good approximation 
for the treatment of pairing correlations is provided by the BCS formalism. Following the 
prescription from Ref. j^j, we use a 5-interaction in the pairing channel, supplemented with 
a smooth cut-off determined by the Fermi function of single-particle energies j^: 

A = ^ _^ ^ie^:-\r-AEr)/f,r ' 

where is the Fermi energy for neutrons (r = n) or protons {r = p). /| is used in the 
evaluation of the pairing density 

Kr{r) = -2^fkUkVkipk{r)^'ipk{r), (7) 

fc>0 

where the summation runs for r = n{p) over neutron (proton) single-particle states. The 
cut-off parameters AEg and fig = AEg/10 are adjusted to the density of single-particle 
levels in the vicinity of the Fermi energy. In particular, the sum of the cut-off weights 
approximately includes one additional shell of single-particle states above the Fermi level 

5^2/fc = iV^ + 1.65iV2/3^ (8) 

A:>0 

where Nj. denotes number of neutrons (protons) in a specific nucleus. The pairing contribu- 
tion to the total energy is given by 

K 
4 

where Ki(p) denotes the strength parameter of the pairing interaction for neutrons (protons). 
Finally, the expression for the total energy reads 

Etot = J [SRMFir) + S;,„{r) + S;,^rir)] dr , (10) 



Elair = / £pairir)dr = \ K\{r)Kr{r)dr , (9) 



and the center-of-mass correction is included by adding the expectation value 

to the total energy, where Pcm denotes the total momentum of the nucleus. 

The principal disadvantage of describing pairing correlations in the BCS approximation is 
that the resulting wave function is not an eigenstate of the particle number operator. More 
precisely, the BCS ground state contains admixtures of particle-number eigenstates with a 
relative spread of order 1 / V^, where N denotes the average number of valence particles. The 
ideal solution, of course, is to perform particle number projection from the BCS state before 
variation. This procedure is technically rather complicated and very much time consuming, 
and therefore it is usual to ernp 
particle number projection [g], Q, 

of local density functionals of the effective interaction, as developed in Refs. 

The Lipkin-Nogami equations are obtained from the variation of the functional 

T=n,p 



m the Lipkin-Nogami (LN) approximation to the exact 
Sj. In this work the LN method is implemented in terms 
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121. 



(12) 



with respect to the single-particle states ipk and the occupation amplitudes Vk- £tot is the 
total energy functional of Eq. flTUl) . The resulting expression for the occupation probabilities 
can be cast into the standard BCS formula 
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A. 



(13) 



where = + 4,\2^tv\ denotes the renormalized single-particle energy, and A,- = Ai^,- -|- 
4A2.T(^r + 1) is the generalized Fermi energy. Ai^,- is determined by a particle number 
subsidiary condition such that the expectation value of the particle number operator equals 
the given number of nucleons. The state-dependent single-particle gaps are defined as the 
matrix elements 

A, = y dr^l{r)^r{r)i,k{r) , (14) 

of the local pair potential 



^r{r) = -VrKr{r) {t = Tl, p) . 



(15) 
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While the quantities Ai,,- represent the Lagrange multiphers used to constrain the average 
particle numbers, the value of the parameters A2,t are determined from 

{HAND 



(16) 



where N2.T is the term of the particle number operator which projects onto two-quasiparticle 



states 



and AiV|^ 



N2,T = 2 ^ UkVk{alal + a-^a^) 



(17) 



A:>0 



A^l^ — (iV|^) denotes its variance. The evaluation of the parameters A2,t is 



described in the Appendix. 

After approximate particle-number projection the total binding energy reads 



(18) 



T=n,p 



The strength of the pairing interaction can be determined by comparing the average pairing 
gaps with empirical gaps obtained from nuclear masses. The average gap is given by the 
summation over occupied states with either the occupation probability f| as the weighting 



factor 



13| 



or the factor UkVk 



li 



(v'A)^ 



{uvA)r 



Efc>o fkvjAk 

Zfc>0 fkUkVkAk 
Z]fc>0 fkUkVk 



(19) 



(20) 



The correspondin g ex pressions for the (approximately) particle-number projected average 



pairing gaps read 
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1^: 



(t;2A)^^ = {v^A)^ + A2,. , (21) 

(m;A)^^ = {uvA), + A2,. . (22) 

The local densities and currents which define the energy functional refer to the intrinsic state, 
and are therefore computed from the Eqs. (E]) - ([5]). On the other hand, the densities used 
to evaluate physical observables, such as the mass quadrupole moment, must correspond to 
the (approximately) particle-number projected state. The LN density is simply com put ed 



by replacing the occupation probabilities vf,, with the LN occupation coefficients 



Wk 



y2y2 y2y2(^y2 _ ^2^2 _ |^^2^2^^2 _ ^2jj2 _|_ 2^2y2 1^^^2^2^2 _ ^4y4j 
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(23) 



where x denotes half of the trace 

x = ^Xk. (24) 

fc>0 

We note that LN-corrected quadrupole moments will be used in the constrained self- 
consistent relativistic mean- field calculations (see Eg. [26i) . 

In this work we only consider even-even nuclei that can be described by axially symmetric 
shapes. In addition to axial symmetry and parity, symmetry with respect to the operator 
^-mjy^ and time-reversal invariance are imposed as self-consistent symmetries. The single- 
nucleon Dirac eigenvalue equation is solved by expanding the large and small components of 
the nucleon spinor ibj in terms of eigenfunctions of an axially symmetric harmonic oscillator 
potential (see Ref. [l?!] for details). 

B. Configuration mixing of mean-field solutions projected on angular momentum 
and particle number 

Correlation effects related to the restoration of broken symmetries and to fluctuations of 
collective coordinates can be taken into account by performing configuration mixing calcu- 
lations of projected states. The generator coordinate method (GCM), which uses a set of 
mean-field states \4>{q)) that depend on a collective coordinate q, provides a very efficient 
procedure for the construction of the trial wave function [18] : 

|v&.) = 5^/a(g,)|0(g,)). (25) 

j 

In this work the basis states \4>{q)) are Slater determinants of single- nucleon states generated 
by solving the constrained relativistic mean-field + LNBCS equations, with the quadrupole 
moment as the generating coordinate q. For an axially deformed nucleus the map of the 
energy surface as a function of deformation is obtained by imposing a constraint on the mass 
quadrupole moment. The method of quadratic constraint uses an unrestricted variation of 
the function 

(/C) + I ((g) - qf , (26) 
where (/C) is the energy functional of Eq. ( |T2|) . (Q) denotes the expectation value of the 
mass quadrupole operator, q is the deformation parameter, and C is the stiffness constant. 

The axially deformed mean-field breaks rotational symmetry, and the particle number 
is only approximately restored with the Lipkin-Nogami procedure, i.e. the basis states 
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10(g)) are not eigenstates of the total angular momentum and particle number operators. 
Therefore, in order to be able to compare model predictions with data, we must construct 
states with good angular momentum and particle number, by performing projections from 
the mean-field plus LNBCS solutions 

l^::^) = E fi^'iiMiKp'p'' \m) ■ (27) 

The particle-number projection operators read 

= — r rf<^„e^(*-^)^" , = — r difpe'^^-^^^^ . (28) 

where N{Z) is the number operator for neutrons (protons), and N{Z) denotes the number 
of neutrons (protons). 

The angular momentum projection operator is defined by 

P'mk = / dQDi:^{Q)R{Q) , (29) 

where the integration is performed over the three Euler angles a, (3, and 7. D'ljj^{Q) = 
e"**'^"c?Mi^(/3)e-*-^'^ is the Wigner function Q, and R{n) = e-'°'^^e-'l^^ye-'^^'' is the rotation 
operator. 

The weight functions fi^{qj) are determined by requiring that the expectation value of 
the energy is stationary with respect to an arbitrary variation 5/^^: 

^}m1,:J (30) 

This leads to the Hill- Wheeler equation 20|]: 

E /"''(^^■) {^'p(<ii)\HP'MKP''p' \m) - K mi)\piiKP''p' m,))) = • (31) 

The restriction to axially symmetric configurations (J^ |0(q')) = 0) simplifies the problem 
considerably, because in this case the integrals over the Euler angles a and 7 can be per- 
formed analytically. In addition, the symmetry with respect to the operator e~™'^y reduces 
the integration interval over the Euler angle (5 from [0,7?] to [0,7r/2]. For an arbitrary 
multipole operator one thus finds 

{m\Qx,P'MKP''P'm,)) = {2J+1) ^^^~^^' sm^,Sko X 

d^r. / d^, / sin/5rf^;o(/?) (0(g,)| gAMe-^^'^-e^(^-^)'^"e^(^-^)'^- |0(g,-)) d(3 .(32) 



We note that this expression is defined only for even values of the angular momentum J. 
The norm overlap kernel 



(2J + 1) 



1 



1)-^ 1 



(27r)2 



2-K p2tt 

d(pn / d(pp X 



tt/2 



sin/3doo(/5) e-'^'ye'(''-''^^"e'^^-^^^- df3 , (33) 



can be evaluated by using the fact, that \(p{qj)) is a product o: 



Slater determinant and the generalized Wick theorem 21 



22 



23 



a neutron- and a proton 



2^: 



n(g„ q,; (3, ^ {<P{q,)\ e'^^'y e'""^^^ |0(g,)) = ±^det Ua,{q^^ qf, /?, V>r) , (34) 
for T = n,p. The overlap matrix is defined as 

J^abiqi, qj; P, ^r) = Ua{qi)Rab{qu qj; P)ubiqj) + Va{qi)Rab{qi, qj; p)vbiqj)e^''^^ , (35) 

where u and v denote the BCS occupation probabilities, and the elements of the matrix R 
read 

Rab{q^, qj; 13) = j ^i(r; g,)e-*^-^-V^b(r; q^)dr . (36) 

We note that the global phase of the overlap in Eq. flMl) is determined by using the procedure 
described in Ref. The details of the evaluation of the matrix R can be found in Ref. 
The Hamiltonian kernel 



= {<P{qi)\HPi,j,P'^P'\ct>{q,)) = 
(2 J + 1) OmqOko 



2tx 



2-K 



■k/2 



in \2 I ^'^^ / ^'^P ^ 

sin/34*(/?) (0(g,)| ^e-^^^-e*(*-^)^"e^(^-^)'^'' |0(g,)) d(3 , (37) 



can be calcu 
densities 21 



ated from the mean-field energy functional Eq. ([T]), provided the modified 
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23|, 



24| 



-r^{r;qi,qj,f3,y^r) = ^Va{qi)vb{qj)e^''^^Aff,^^{qi,qj; (3) 



X 



Mf;qi){-h^ + m)e '^■^yipb{r;qj) , 

J2 Mqi>b{qjy^^Ka\qh qr, q^e-'^^yMr; qj] 



(38) 
(39) 



f^ir;qi,qj,l3,ipr) = ^Va{qi)vb{q,)e^'''-M^^\qi,qf,P)iJa{r-qi)-i^e-''^^y^b{r]qj), (40) 



a,b 



are used when evaluating the expression 



The computational task of evaluating the Hamiltonian and norm overlap kernels can be 
reduced significantly if one realizes that states with very small occupation probabilities give 
negligible contributions to the kernels, and hence such states can be excluded from the 



calculation 



23 



24|. 



For an even number of particles, the integration interval in Eq. fl28|) can be reduced to 
[0, tt] using the symmetries of the i nteg rand. Furthermore, the integrals can be discretized 
by using the Fomenko's expression [25| 



'^ = lX^e^(^-^)-", ¥.„ = ^n, (42) 

n=l 

with L points in the expansion. In order to avoid numerical instabilities which might arise 
at V9 = I when the occupation probability of a state is exactly 0.5, an odd number of points 
must be used in the Fomenko's expansion |26:|. The Gauss-Legendre quadrature is used for 



the integration over the Euler angle /3. We have verified that already L = 9 points both 
for protons and neutrons in the Fomenko's expansion, and 13 points in the integral over (3, 
produce numerically stable results. 
The Hill- Wheeler equation ( l3Ti) 



presents a generalized eigenvalue problem. Thus the weight functions fa{qi) are not orthog- 
onal and cannot be interpreted as collective wave functions for the variable q. The standard 



27l | is to re-express Eq. fH3|) in terms of another set of functions, gaili), defined 



procedure 
by 

j 

With this transformation the Hill- Wheeler equation defines an ordinary eigenvalue problem 

J2T^'iq^,qj)9iiq,) = E^giiq,), (45) 
j 

with 

n■\q^,qJ) = E(-^')"'^'(^-^'^)^'(^'^'*)(-^')"'^'(^''^^) • (46) 



k,l 
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The functions QaiQi) are orthonormal and play the role of collective wave functions. For a 
more detailed description of this particular implementation of the Hill- Wheeler equation, we 
refer the reader to Ref. 

For completeness we also include the expressions for physical observables, such as tran- 
sition probabilities and spectroscopic quadrupole moments |2S]. The reduced transition 
probability for a transition between an initial state (Jj, ai), and a final state {Jf, af), reads 

2 



B{E2-Jiai Jfaf] 



e2 



qf,qi 



2J^ + l 

and the spectroscopic quadrupole moment for a state (Ja) is defined 



(47) 



Q'''V,c^) = eJ^ (J'^ J {Jq^\ \Q2\ \Jq,) fM ■ (48) 

\ ^ ^ ^ / <1^A3 

Since these quantities are calculated in full configuration space, there is no need to introduce 
effective charges, hence e denotes the bare value of the proton charge. In order to evaluate 
transition probabilities and spectroscopic quadrupole moments, we will also need the reduced 
matrix element of the quadrupole operator 



{JJqJ\m\J^q^) = (2J, + 1)(2J; + 1) V [ ^ J/ \ 1 

^ \-fi fi / (2^) Jo Jo 



'•271 r2TT 



tt/2 







smpd'j;,{P) (0(g^)|g2Me-^^"^e^(^-^)'^"e^(^-^)^'' |0(g,)) f//3.(49) 



III. ILLUSTRATIVE CALCULATIONS 

The intrinsic wave functions that will be used in configuration mixing calculations are ob- 
tained as solutions of the self-consistent RMF+LNBCS equations, subject to constraint on 
the mass quadrupole moment. As in the first part of this analysis we use the relativistic 
point-coupling interaction PC-Fl ^ in the particle-hole channel, and a density-independent 
5-force is the effective interaction in the particle-particle channel. The parameters of the 
PC-Fl interaction and the pairing strength constants V„ and Vp have been adjusted simul- 
taneously to the nuclear matter equation of state, and to ground-state observables (binding 
energies, charge and diffraction radii, surface thickness and pairing gaps) of spherical nuclei 
{4], with pairing correlations treated in the BCS approximation. However, since the present 

11 



analysis includes the LN approximate particle number projection, the pairing strength pa- 
rameters have to be readjusted. By comparing the projected average pairing gaps {uvA)i^^^ 
and the BCS pairing gaps {uvA)i^'"^\ we find that the neutron pairing strength should be 
reduced from Vn = —308 MeV to Vn = —285 MeV, and the proton pairing strength from 
Vp = —321 MeV to Vp = —260 MeV. The average pairing gaps for two isotopic and two 
isotonic chains are shown in Fig. [H We notice that by readjusting the strength parameters 
Vn and Vp, a good agreement between projected average pairing gaps and the BCS average 
pairing gaps is obtained except, of course, for magic numbers of neutrons or protons, for 
which an unphysical collapse of pairing correlations is found in the BCS approximation. 

In order to illustrate the importance of including particle-number projection in the de- 
scription of specific structure effects, we will compare the results of GCM calculations with 
those obtained using the model that was developed in the first part of this work [1], and 
in which the intrinsic wave functions are generated from the solutions of the constrained 
relativistic mean-field + BCS equations, without the Lipkin-Nogami approximate particle- 
number projection. In that case the correct mean-values of the nucleon numbers are restored 
by modifying the Hill-Wheller equations to include linear constraints on the number of pro- 
tons and neutrons (see Eqs. (54) and (55) of Ref. [l|). Therefore, in the following subsections 
AMP will indicate that only angular momentum projection has been carried out before GCM 
configuration mixing (i.e. the model of Ref. |l| is used), and PN&AMP will denote the re- 
sults of GCM calculations which include both the restoration of the particle number and 
rotational symmetry. 

In this section illustrative configuration mixing calculations are presented for ^^Mg, ^^S 
and ^^Ar. We choose these nuclei because the results can be directly compared with extensive 
GCM studies performed using the nonrelativistic Skyrme and Gogny effective interactions. 



In the analyses in which the nonrelativistic zero-range Skyrme interaction was used [23|, |29| , 
both particle number and angular momentum projections were performed. The simultane- 
ous projection on particle number and angular momentum is computationally much more 

ore only angular momentum 



28 



30 



3l| . Both approaches 



demanding in the case of a finite-range interaction, and there 
projection was performed in the studies with the Gogny force 
are interesting for the present analysis, because by comparing the results one can deduce 
which effects can be attributed to particle number projection, and which originate from the 
differences in the properties of the effective interactions. 
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The constrained RMF equations are solved by expanding the Dirac single-nucleon spinors 
in terms of eigenfunctions of an axially symmetric harmonic oscillator potential. In order 
to keep the basis closed under rotations, the two oscillator length parameters b± and bz 



have always identical values 



32 



331]. In addition, to avoid the completeness problem in 



subsequent configuration mixing ca^ 



culations, the same oscillator length is used for all values 



of the quadrupole deformations 33||. Since we consider relatively light nuclei in this work. 



it is sufficient to use ten oscillator shells in the expansion (see Ref. l|] for details). 



A. 24Mg 

The low-energy spectrum of ^^Mg displays a typical rotational structure, with data on the 



ground-state band extending up to angular momentum I = 8H 



34 



35| . The spectroscopic 



quadrupole moment of the 2f state: —16.6(6) e fm^ |35|, indicates a large prolate defor- 
mation. The 0^ state is found at high excitation energy (6.432 MeV), and this means that 
shape-coexistence effects should only play a minor role in the description of the ground-state 



band. Properties of Mg have been studied with the GCM using the Skyrme SLy4 



and the Gogny DIS 



28 



23|, 



effective interactions. 

In Fig. [2] we display the pairing energy (upper panel), and the total RMF binding energy 
curve (lower panel) of ^^Mg, as functions of the mass quadrupole moment. The Lipkin- 
Nogami procedure has not been implemented at this stage, and consequently pairing cor- 
relations vanish in a broad region of deformations around the deformed first minimum of 
the potential energy surface. Since the moment of inertia of a rotational band is reduced 
in the presence of pairing correlations, dynamical pairing effects could be important in the 
description of the ground-state band of ^"^Mg. 

The GCM excitation energies and the resulting transition probabilities for the ground- 
state band, calculated with the PC-Fl effective interaction, are shown in Fig. [3l The results 
of the AMP and PN&AMP configuration mixing calculations are compared with the data. 
As expected, the inclusion of dynamical pairing effects reduces the moment of inertia, but 
the resulting spectrum is much too spread out compared to experiment. This is a well 
known problem, related to the fact that we project particle number and angular momen- 
tum only after variation, rather than performing the projections before variation [36]. It 
has been shown that in the latter case rotational bands with larger moments of inertia are 
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obtained [37|, |38j, provided that the model geometry allows for the alignment of nucleon 
angular momenta. Since the full projection before variation is technically and computation- 
ally much more complex, it has been seldom used in realistic calculations. We note that 
one possible improvement of the present model would be to project, for each value of the 
angular momentum J, the cranked mean-field wave functionswhich, in addition to the mass 

391. This extension has not 



quadrupole moment, are also constrained to have {J^) = J [3, 
been included in the present analysis. 

The transition probabilities, as well as the calculated spectroscopic quadrupole moment 
Qspeci'^i) = —16.56 e fm^, are in very good agreement with the data. The AMP GCM 
results are very similar to those obtained by using the Gogny DIS force (see Figs. 11 and 



13 in Ref. 



281]), whereas the PN&AMP spectrum is close to the one calculated with the 



Skyrme SLy4 interaction (see Fig. 6 in Ref. 23l|). 

The amplitudes of the PN&AMP GCM collective wave functions \gi\'^ of the ground-state 
band, are plotted in Fig. H] as functions of the mass quadrupole moment. It is interesting 
to note that only the 0^ state contains significant admixtures of oblate deformed shapes, 
whereas the amplitudes of states with J > 2 are concentrated in the prolate well. The same 
result has also been obtained with the Skyrme SLy4 interaction (see Fig. 4 of Ref. [23]). 



B. 



32(1 



In recent years a number of theoretical and experimental studies of the structure of ^^S 
have been reported. This nucleus is particularly interesting because the excitation energies 
of the low- lying 0^, 2'^ and states correspond to those of a typical spherical vibrator, 
whereas on the other hand, the quadrupole moment of the first 2+ state is large and negative, 
indicating large dynamical prolate deformation [4^. ^^S is among the best studied nuclei in 
the sd shell, and both the energies and lifetimes of many states up to ~ 10 MeV excitation 



energy are known 
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42|. 



Several modern theoretical approaches have been used in recent studies of normally de- 
formed (ND), and superdeformed (SD) configurations in ^^S: the shell model with the univer- 
sal sd-shell Hamiltonian 43|], the semimicroscopic algebraic cluster model 4J], and various 
extensions of the self-consistent mean-field framework. They include the cranked Hartree- 



Fock 
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471], and Hartree-Fock-Bogoliubov method [30] for the description of the SD 
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configuration, and the generator coordinate method with the Skyrme SLy6 |29| . and Gogny 
DIS effective interactions, in the analysis of both ND and SD configurations^, 13 ■ 

In Fig. O the pairing energy (upper panel), and the total RMF binding energy curve 
(lower panel) of ^^S, are plotted as functions of the mass quadrupole moment. The effective 
interaction is PC-Fl, and pairing correlations are treated in the BCS approximation, without 
particle number projection. Consequently, we notice the collapse of the pairing energy both 
in the ND and the SD minima. This means that restoring the particle number symmetry 
could have an important effect in the description of both ground-state and SD bands. 

The results of the PN&AMP and AMP GCM calculations for ^^S are shown in the left 
and right panels of Fig. El respectively. In the left (right) panel we also include the mean- 
field MF+LNBCS (MF+BCS) binding energy curve (dotted curves). In both cases the 
MF energy curve displays a spherical minimum, and an additional shallow minimum at 
large deformation g ^ 4b with excitation energy ^ 11 MeV. The occurrence of almost 
degenerate oblate and prolate minima in the J = 0^ projected energy curve, symmetrical 
with respect to the spherical configuration, is a feature common to all nuclei for which the 



mean-field calculation predicts a spherical ground state (see, for instance, Refs. (29|, |48|). 

The superdeformed minimum is more pronounced in the calculation without particle 
number projection, both for the MF+BCS curve and for the angular momentum projected 
energy curves. The GCM superdeformed band is calculated at somewhat lower excitation 
energies in the AMP case, and this is because in the PN&AMP calculation pairing corre- 
lations do not vanish in the SD minimum. The energies of the GCM states are plotted as 
functions of the average quadrupole moment 

Q = ilk) = ^gl{qj)qj ■ (50) 
j 

The overall structure of energy levels, and in particular the two-phonon triplet 0^, 2^, 4^, 
is much better described when particle number projection is included (see Tables |T] and |Tl] 
for a comparison with experimental levels). Similar results for the spectra of '^^S have also 



been obtained in the GCM analyses of Refs. 29|, ISOj. The calculation with the zero-range 



SLy6 effective interaction, including particle number projection, reproduces the structure of 
the two-phonon triplet [29] . On the other hand, the results for the ND states obtained with 



the finite-range Gogny force, but with only angular momentum projection 30|, are similar 



to those shown in the right panel of Fig. [6l i.e. an additional low-lying 0"*" state is predicted 

15 



by the calculation. 

In Tab. [J we list the PN&AMP GCM excitation energies, the spectroscopic quadrupole 
moments, and the E2 (J — > J -2) transition probabilities for the 2^ state, and for the two- 



4 



and the available 



phonon triplet 0^, 2^, 4^, in comparison with the results of Ref. 
data. The results obtained with both effective interactions, the nonrelativistic SLy6 and 
the relativistic PC-Fl, are in qualitative agreement with the data. The excitation energies 
of the two-phonon triplet are calculated at approximately twice the energy of the 2f state. 
We note, however, that the predicted quadrupole moments for the 2^ state are small and 
positive, whereas the experimental value points to a much larger and prolate deformation 
of this state. 

A small prolate deformation for the ground-state band is only obtained in calculations 
which do not include particle number projection (right panel of Fig. [6]). In Tab. [11] we 
compare the AMP GCM excitation energies and B(E2, J — > J -2) values obtained with PC- 



Fl, with the corresponding quantities calculated with the Gogny DIS interaction [SOj. The 
results are very similar, with the largest differences in the B(E2) values for the transitions 
2f 0^ and 2^ O3 . This can be attributed to the smaller values of the quadrupole 
moments of the 2^ states predicted by the PC-Fl interaction. For instance, the spectroscopic 
quadrupole moment of the 2f state calculated with the Gogny force —13.29 e fm^ is close 
to the experimental value: —14.9 e fm^ [s^, whereas the one predicted by the PC-Fl 
interaction: —3.5 e fm^, is much too small. 

Although the present version of the RMF plus GCM model is not optimal for the study of 
moments of inertia of superdeformed rotational bands, because it does not include cranked 
wave functions, nevertheless it can be used to investigate the stability of the SD intrinsic 
configurations against quadrupole fluctuations at low angular momentum. In Fig. [6] it is 
shown that both the PN&AMP and AMP GCM calculations predict an SD band at large 
deformation q ^ 46. In the left panel of Fig. [7] we plot the energy differences AE{J) = 
E{J) — E{J — 2), as functions of the angular momentum of the SD band. Since without 
particle number projection pairing correlations vanish in the SD minimum, the AMP GCM 
calculation overestimates the moment of inertia of the SD band, and the resulting spectrum 
is more compressed. The moment of inertia is reduced with the inclusion of dynamical 
pairing effects. In the right panel of Fig. [7] we include the corresponding E2 transition 
probabilities. Both the relative excitation energies and the B(E2) values for the SD band 
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are consistent with the results reported in Ref. [29| (Tab. VIII), and Ref. [30| (Fig. 3). We 
note, however, differences in the predicted position of the SD band-head. In our PN&AMP 
GCM calculation the SD band-head is found at 8.9 MeV, almost 3 MeV lower than the 
value calculated with the SLy6 interaction in Ref. j^j. The AMP GCM calculation predicts 
the head of the SD band at 7.6 MeV, comparable to the value obtained with the Gogny 
DIS interaction in Ref. 30(] (8.9 MeV). A possible reason for the large difference between 



the PC-Fl and Gogny DIS interactions on one hand, and the SLy6 Skyrme force on the 
other, can be identified already at the mean-field level. While the energy curves obtained 
with the PC-Fl and Gogny interactions (see also Fig. 2 in Ref. ^]) exhibit a second, 
superdeformed shallow minimum at g ~ 46, only a shoulder in the potential energy curve at 
super deformation is predicted by the SLy6 interaction (see Fig. 13 in Ref. |29l|). 



C. 



36 



Ar 



In the last example we present the GCM analysis of the low angular momentum structure 
of 36 Ar, one of the lightest nuclei in which a superdeformed structure has been studied 
experimentally. The excitation energies and E2 transition probabilities for the prolate SD 



band have been recently measured up to the angular momentum / = 16H 49 . The properties 
of the oblate ground-state band were determined over a decade ago 34|. A number of 
theoretical analyses of the structure of ^^Ar include the cranked Nilsson-Strutinsky and the 



shell model 



49 



Bogoliubov model 



50|, the projected shell model |51|, the self-consistent cranked Hartree-Fock 



3l| . and the generator coordinate method with Skyrme SLy6 [29| and 



Gogny DIS [31] interactions. 

The results of the PN&AMP and the AMP GCM calculations are shown in the left and 
right panels of Fig. [HI respectively. In both cases the mean-field binding energy curve is also 
included (dotted curves). The BCS mean-field energy curve displays an oblate minimum, 
rather flat in the region of deformation: — 1 b < q < 0.5 b. In addition, a shallow minimum is 
found at larger deformation q ^ 2.8 b, at -E^, ~ 9 MeV above the ground-state minimum, and 
a shoulder is predicted at still larger deformation g ^ 5 b. A similar mean-field potential 
energy surface is obtained with the Gogny DIS interaction (see Fig. 1 in Ref. Ml]), but 
the SD minimum is calculated ~ 1 MeV lower than with PC-Fl. We note that in the BCS 
approximation pairing correlations vanish both in the ND and SD minimum, and this means 



17 



that the restoration of particle number can produce sizable effects in the ground-state and 
SD bands. The mean-field energy curve which includes the approximate particle number 
projection (LNBCS), displays only a weak shoulder, rather than a minimum at deformation 
g ^ 3 b. This curve is in qualitative agreement with the LNBCS binding energy curve 
calculated with the Skyrme SLy6 effective interaction (see Fig. 10 in Ref. ^]), although 
in the latter case the shoulder occurs at smaller deformation g ^ 1.8 b, and much smaller 
excitation energy 5 MeV. Both in the PN&AMP and the AMP calculations, the J = 
angular momentum projected energy curves display two well developed low-lying minima 
on the oblate and prolate side. The oblate minimum corresponds to the ground state. An 
additional SD minimum occurs on the AMP J = energy curve, whereas the corresponding 
PN&AMP curve displays only a plateau. The angular momentum projected energy curves 
with J > 2 exhibit well developed SD minima both in the AMP and PN&AMP calculations. 
Finally, the energies of the resulting GCM states for angular momenta J < 8 are included 
as functions of the average quadrupole moment, defined in Eq. ( !50l) . The low-spin GCM 
spectra contain an oblate ND ground state band, and a prolate SD band. 

In Tab. IIIII we compare the PN&AMP GCM excitation energies, the spectroscopic 
quadrupole moments, and the B(E2, J —>■ J-2) values for the ground-state band, with the 
results obtained with the Skyrme SLy6 interaction in Ref. I29I, and with available data. 
The theoretical results are in qualitative agreement. Although the energy spectrum is de- 
scribed somewhat better by the PC-Fl interaction, both calculations clearly underestimate 
the moment of inertia of the ND ground-state band. In Tab. IIVI the AMP GCM excitation 
energies and the B(E2, J J-2) values for the ND band, are shown in comparison with 
those calculated with the Gogny DIS interaction in Ref. 31|- The results obtained with the 
PC-Fl and Gogny effective interactions are very similar, and a considerable increase of the 
moment of inertia as compared to the PN&AMP GCM spectra, is caused by the collapse of 
pairing correlations in the ND minimum. 

In the left panel of Fig. Owe plot the energy differences AE{J) = E{J) — E{J — 2), and 
in the right panel the BE2 values, as functions of the angular momentum of the SD band, 
in comparison with data. Cranked Hartree-Fock-Bogoliubov calculations performed with 
the Gogny interaction [sijl, have shown that rather strong triaxiality effects appear already 
at zero-spin and, as a result, predicted a steady decrease of deformation with increasing 
angular momentum. Being restricted to axially symmetric shapes, the model that we have 
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developed in this work cannot take these effects into account. This is clearly reffected in 
the pronounced discrepancy between the calculated and experimental B(E2) values. Similar 
results have also been obtained with the SLy6 (see Tab. V in Ref. 129|), and Gogny effective 
in.e.ac«o„. (.e Fi. 3 in Ref. E). Ju. U.e . t.e ca.e of ... we W con.de.ed in 
the previous section, the transition energies and the B(E2) values do not crucially depend 
on the effective interaction in the particle-hole channel. On the other hand, the predicted 
positions of the SD band- head differ significantly: for the Gogny DIS interaction the band- 
head is at 7.5 MeV, whereas the SLy6 interaction predicts a lower value of 5.9 MeV. The 
positions of SD band-head obtained in the PN&AMP and AMP GCM calculations with the 
PC-Fl interaction are 9.2 MeV and 9.4 MeV, respectively. 

Finally, the amplitudes of the collective PN&AMP GCM wave functions \gi\^ for the ND 
(0^, 2|, 4^ and 6^), and the SD (O;]", 2 J, 4^ and and 8^^) rotational bands in ^^Ar are 
plotted in Fig. [TOl We notice that, except for the ground state 0^, the amplitudes of the 
states of the ND band are principally concentrated in the oblate minimum, whereas in the 
panel on the right the amplitudes of the wave functions are strongly peaked in the prolate 
superdeformed minimum. 



IV. SUMMARY AND OUTLOOK 



In Re 
theory {2 



'. 11] and in this work we have extended the very successful relativistic mean-field 
3|] to explicitly include correlations related to the restoration of broken symmetries 
and to fluctuations of collective coordinates. We have developed a model that uses the 
generator coordinate method to perform configuration mixing of angular-momentum and 
particle-number projected relativistic wave functions. The geometry is restricted to axially 
symmetric shapes, and the intrinsic wave functions are generated from the solutions of 
the relativistic mean-field + Lipkin-Nogami BCS equations, with a constraint on the mass 
quadrupole moment. The single-nucleon Dirac eigenvalue equation is solved by expanding 
the large and small components of the nucleon spinor in terms of eigenfunctions of an axially 
symmetric harmonic oscillator potential. The current implementation of the model employs 
a relativistic point-coupling (contact) nucleon-nucleon effective interaction in the particle- 
hole channel, and a density-independent 5-interaction in the particle-particle channel. 
We have performed several illustrative calculations which test our implementation of the 
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GCM with simultaneous particle-number and angular-momentum projection. The results 
have been compared with those obtained employing the relativistic GCM model developed 
in the first part of this work [l| , which does not include particle-number projection, but only 
conserves the number of particles on the average, and with results that were obtained using 
the corresponding non-relativistic models based on Skyrme and Gogny effective interactions. 
In this work the low- lying spectra of ^^Mg, ^^S and ^^Ar have been analyzed. The results 
of GCM configuration mixing calculations both without and with particle-number projec- 
tion, have been compared with those obtained with the angular-momentum projected GCM 
based on the Gogny DIS effective interaction, and with the particle-number and angular- 
momentum projected GCM based on the Skyrme SLy4 and SLy6 effective interactions. The 
principal features of the spectra calculated without (Gogny interaction) and with particle- 
number projection (Skyrme forces), are very well reproduced in our GCM calculation with 
the PC-Fl relativistic effective interaction. These include the deformations and moments 
of inertia of yrast bands, and the occurrence and structure of superdeformed bands. Some 
differences between the predictions of the three models, as for instance the position of the 
head 0"*" of the superdeformed band, can be attributed to the gaps in the mean-field single- 
nucleon spectra calculated with the different effective interactions. We have also shown that 
dynamical pairing effects play an important role in the description of the low-energy spectra. 
In particular, we find a pronounced effect of particle-number projection on the moments of 
inertia of the ground-state and superdeformed rotational bands. Another example is the 
two-phonon triplet 0^, 2^ and in ^^S, which can be reproduced by GCM calculations 
only when dynamical pairing effects are taken into account. 

There are, of course, many possible improvements and extensions of the present imple- 
mentation of the relativistic GCM model. Perhaps the most obvious is the extension to 
shapes that are not constrained by axial symmetry. The inclusion of triaxial deformations 
is in principle straightforward but, because it requires an enormous increase of computa- 
tional capabilities, not feasible at present. The second major problem is that our GCM 
configuration mixing calculations correspond to a projection after variation. A more general 
variation after projection is far too complicated to be used in realistic calculations at the 
present stage. A possible improvement, however, is to generate the GCM basis functions, for 
each value of the angular momentum, by performing cranking RMF-I-LNBCS calculations 
with the additional constraint {Jx) = a/ J{J + !)• This would automatically increase the 
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moments of inertia of rotational bands, and therefore produce spectra in better agreement 
with experiment. The extension to odd-A nuclei necessitates the breaking of time- reversal 
invariance in the wave functions and, therefore, the explicit inclusion of currents in the 

energy-density functional, and the corresponding time-odd fields in the single-nucleon Dirac 

i 

equation. Finally, let us emphasize again one the conclusions of Ref . , namely that those 
correlations which are explicitly treated in the GCM configuration mixing, should not be 
contained in the effective interaction in an implicit way, i.e. by adjusting the parameters 
of the interaction to data which already include these correlations. Therefore, before the 
present version of the relativistic self- consistent GCM is applied in realistic calculations, we 
need to adjust a new global effective interaction which will not contain symmetry breaking 
corrections and quadrupole fluctuation correlations. 



APPENDIX: EVALUATION OF THE PARAMETER A2 IN THE LIPKIN- 
NOGAMI METHOD 



The variation of the auxihary functional /C Eq. ([T2|) can be expressed as 







(A.l) 



for the ground-state expectation value, with N2 as the two-quasiparticle part of the particle- 
number operator Eq. (fT7|) . Since we do not consider proton- neutron pairing, the equations 
for protons and neutrons separate, and the isospin index q can be omitted in the following 
derivation. If one introduces the shifted many-body state 



10 



|0) 



|o) = 10 



(A.2) 



the following expression for A2 is obtained from the variational condition Eq. (lA.ip 111 ] 

dU^\H\i) 



X2 



€=0 



dUo\m\0 



(A.3) 



The denominator of this equation is evaluated using Wick's theorem 



dUo\N'\o 



32 



(A.4) 
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The principal advantage of this particular formulation of the Lipkin-Nogami scheme is that 
one can relate it to the theoretical framework of energy density functionals 



(0| H 10 ^ £^^^ = £[p^^\ k^^\ k<^^] 



with the following definition of the shifted densities 



pi\^ = maW\0, 4? = (0| arable) 



(A.5) 



(A.6) 



The second derivative of the energy functional reads 
S£ „n . S£ ^r, ^ S£ 



= Tr 
+ 2Tr Tr 



where the trace implies integration and summation over all coordinates. The term with a 



single trace vanishes in the BCS ground state [ll| and, since we use volume pairing in the 
energy functional (see Eq. [9]), the mixed derivative terms which contain p and k, or p and 
k*, also vanish. With the definition of the response densities p, k and k* 



p = —id^p^^^ , K = —idfk^^^ 
one finally obtains 



and K* = —idfk''^'** 



5=0 



5=0 



= -Tr Tr 



5^£ . . , 5-'£ 

-piP2 + 2 



Spi6p2 6k6k* 
The RMF energy functional considered in this work consists of three terms: 
i) the kinetic energy 



(A.8) 



(A.9) 



£..^j:Jclr.lMr)H.V-,m)Mr) 



(A.IO) 



ii) the field energy 

+ -i^3tv{3tv)i, + —3tv/^{3tv)i. + ^Pts + —pTS^pTsj , (A.ll) 



and iii) the Coulomb interaction term 

^2 



£ 



c 




drdr 



,Ppir)pp{r') 



r — r 



l\2 



(A.12) 
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The response densities which appear in Eq. f lA.QP are given by 

pSW = 5^«XV^.(r)^,(r) (A.13) 

k 

•^mW = Y.'^kvlMrHMr) , (A.14) 

k 

where the summation runs for r = n{p) over neutron (proton) single-particle states. The 
functional derivative of S field reads 

disfieid = / [("^ + 2/^sP5 + 375P| + aTs)pVpV + (Ss + 6ts)pVApV 

+ K + 375JmJ^ + «Tv)j:?" + (5y + 5Ty)j:A/1 , (A.15) 

where ps and denote the scalar density and baryon current, respectively , and and 
are the corresponding neutron (proton) densities and currents. For protons there is an 
additional contribution from the Coulomb interaction: 




To evaluate the contribution of the pairing energy Eq. ([9]) to the second derivative of the 
energy functional, one needs the non-hermitian response pairing tensor 



k^{r) = -4j2fk4vk^Pi{r)Mr) (A.17) 

fc>0 

k^*ir) = AY,fk^kvl4ir)Mr) ■ (A.18) 



fc>0 

This leads to a rather simple expression 



dlsS-r = Y I drk-K-*. (A. 19) 



Finally, by inserting Eqs. flA.lSp . flA.16p . and flA.lQP into Eq. flA.3p . one obtains the Lipkin- 
Nogami parameter A2,r- 
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FIG. 1: Comparison of the average pairing gaps calculated with a simple BCS, and the Lipkin- 
Nogami plus BCS approximations. For the Z = 50 and Z = 28 isotopes the neutron pairing gaps 
are shown, whereas for the N = 50 and N = 82 isotones the proton average pairing gaps are 
compared. 
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FIG. 3: The ground-state rotational band in ^^Mg. The B(E2) values are in units of e^fm^. The 
GCM spectra calculated with (PN&AMP), and without (AMP) particle number projection, are 
compared with the experimental ground-state band 34 1 
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FIG. 6: The energies and the average quadrupole moments of the GCM states in ^^S, plotted 
together with the corresponding angular momentum projected energy curves. GCM calculations 
with (left panel), and without (right panel) particle number projection are compared. The mean- 
field binding energy curves are also included in the figure (dotted curves). In both panels zero 
energy corresponds to the minimum of the J = projected energy curve. 
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FIG. 7: Transition energies AE{J) = E{J) - E{J - 2) (left panel), and the B(E2, J-^ J-2) values 
(right panel), for the SD configuration in ^^S, as functions of the angular momentum. 
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FIG. 8: The energies and the average quadrupole moments of the GCM states in ^^Ar, plotted 
together with the corresponding angular momentum projected energy curves. GCM calculations 
with (left panel), and without (right panel) particle number projection are compared. The mean- 
field binding energy curves are also included in the figure (dotted curves). In both panels zero 
energy corresponds to the minimum of the J = projected energy curve. 
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FIG. 9: The energy differences AE{J) = E{J) - E{J - 2) (left panel), and the B(E2, J-2) 
values (right panel), for the SD band in ^^Ar, as functions of the angular momentum. Results of 
GCM calculations with and without particle number projection are compared with the available 
data 49 
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FIG. 10: The amplitudes of the PN&AMP GCM eohective wave functions \g(\'^ for the ND ground- 
state band (0]^, 2^ , A'l and 6^), and the SD band (0^, 2^, 4^ and Q'^ and 8]^) in ^^Ar, as functions 
of the mass quadrupole moment. 
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TABLE I: Excitation energies, spectroscopic quadrupole moments and B(E2, J J-2) values for 
the one-phonon (2^) state, and the two-phonon triplet (0^, 2^ and 4^'") in '^^S. GCM calculations 



performed with particle number projection are compared with the results reported in Ref. 
and with available data. 



state 


(MeV) 


Q specie frn^) 


Transition 


BE2{e'^fm'^) 




This work Ref. [29] Exp. 


This work Ref. [29] Exp. 




This work Ref. [29] Exp. 




0.0 0.0 0.0 










2.81 3.22 2.23 


5.8 2.3 -14.9 


2+^0+ 


94 38 61 


0+ 


6.33 6.32 3.78 




0+ ^ 2+ 

U2 


37 144 72 




6.33 7.04 4.28 


-3.0 -0.7 


2+ 2+ 
2+^0+ 
2+^0+ 


131 157 54 
0.134 0.02 11 
10.7 2.8 




6.61 7.35 4.46 


-1.2 11.7 


4+ 2+ 


140.5 94 72 



TABLE IL Excitation energies and B(E2, J J-2) values for the two lowest bands in ^^S. The 
results of GCM calculations performed without particle number projection are shown in comparison 
with those of Ref. 30(] . 



state 



2^ 
6^ 



03^ 



4 



(MeV) 



This work Ref. [30] 



0.0 
2.181 
5.395 
9.661 



0.0 
2.107 
5.825 
10.962 



3.24 3.778 
4.832 4.282 
9.213 9.097 



Transition 



BE2{e^fnv' 



This work Ref. [30 



66.2 72.3 
102.9 119.8 
146.1 142.8 



33.3 
121.1 



58.0 
132.2 
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TABLE III: Excitation energies, spectroscopic quadrupole moments, and B(E2, J J-2) values 
for the ground-state band in ^^Ar. GCM calculations performed with particle number projection 
are compared with those of Ref . {29^ , and with available data. 



state 


(MeV) 


Q specie fm^) 


Transition 






This work Ref. [29] Exp. 


This work Ref. [29] Exp. 




This work Ref. [29] Exp. 


0^ 


0.0 0.0 0.0 








2^ 


2.26 2.8 1.97 


13 13 


2+^0+ 


79 44 60± 6 


4^ 


6.40 7.43 4.41 


14 12 


4+^2+ 


129 103 


6+ 


13.94 13.65 9.18 


9.9 -1.3 


6+^4+ 


152 93 



TABLE IV: Excitation energies and B(E2, J — > J-2) values for the ground-state band in ^^Ar. 
GCM calculations performed without particle number projection are compared with the results 
reported in Ref. [311 ]. 



state 




(MeV) 


Transition 


BE2{e^fm'^) 




This work Ref. [31] 




This work Ref. [31j 


0^ 


0.0 


0.0 








2| 


1.54 


1.45 






74.8 72.1 


4^ 


4.99 


4.54 


4^- 


-2+ 


114.7 102.3 


6^ 


12.15 


10.01 




-4+ 


142.4 112.8 
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